EXTENDING THE RANGE OF ERROR ESTIMATES FOR RADIAL 
APPROXIMATION IN EUCLIDEAN SPACE AND ON SPHERES 



R. A. BROWNLEEtt, E. H. GEORGOULISt , AND J. LEVESLEYt 

Abstract. We adapt Schaback's error doubling trick [R. Schaback. Improved error bounds 
for scattered data interpolation by radial basis functions. Math. Comp., 68(225):201-216, 1999.] 
to give error estimates for radial interpolation of functions with smoothness lying (in some sense) 
between that of the usual native space and the subspace with double the smoothness. We do this for 
both bounded subsets of M.'^ and spheres. As a step on the way to our ultimate goal we also show 
convergence of pseudoderivatives of the interpolation error. 

Key words, multivariate interpolation, radial basis functions, error estimates, smooth functions 

AMS subject classifications. 41A05, 41A15, 41A25, 41A30, 41A63 

1. Introduction. In this paper we are interested in extending the range of ap- 
phcabihty of error estimates for radial basis function interpolation in Euclidean space 
and on spheres. Let be a subset of R'*, or the sphere. Let d{x,y) denote the dis- 
tance between two points in ^l. Let F C be a finite set of points, and measure the 
fill-distance of y in with 

h{Y,n) sup Tain d{x , y) . 

Given a univariate function <f> defined either on R-|_ or [0, tt], depending on whether 
we are in Euclidean space or on the sphere, we form an approximation 

Sl{x) = ^ ay4>{d{x,y)). 

If the coefficients ay, y £Y ^ are determined by the interpolation conditions 

Sliy) = /(y), for y e Y, 

we refer to Sj^ as the 0-spline interpolant to / on y. 

We will be approximating functions / G IK^ , a Hilbert space of functions which 
depends on the function (jj — the so-called native space. Later we will be more explicit 
about this space of functions. With this Hilbert space we have a inner product ( • , ■)^, 
with associated norm || • \\^. We will require the following useful orthogonality and 
consequent Pythagorean property; see, e.g., [TT| I13j. 

Proposition 1.1. Let be the (j)-spline interpolant to f on the point set Y C fl. 
Then, for all f G JC^, 

2- \\f\\l-\\f-sl\\l + \\sj\\l. ^ 

The usual error estimate for 0-spline intcrpolants is of the form 
\fix) - Slix)\ < P{x,Y,mf - S^h, 
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where estimation of P{x, Y, (p) — the so-cahed power function — leads to error estimates 
for interpolation in terms of the fill-distance h{Y,il). For the archetypal function in 
J-Ccf, we can say no more than ||/ — S'^Wif, ^ as h{Y,n) — » 0. However, if / has 
double the smoothness (in some sense to made clear later) than the typical function, 
then Schaback 19J has shown how to double the convergence order of the 0-spline 
interpolant. 

We show how to get improved orders of convergence when the target function, /, 
has less smoothness than Schaback requires, but more smoothness than the typical 
function. We shall be doing this on the sphere (though this can be easily generalised 
to other two-point homogeneous spaces) in ^ and in Sobolev spaces on Euclidean 
space in SJ31 An intermediate result in both cases is to prove approximation orders for 
pseudoderivatives of the interpolant. We will define this notation at the appropriate 
place in each of the following sections. 

In each case we shall be concerned with the practical scenario in which Y consists 
of a finite number of points. Foregoing this assumption is of theoretical interest. In 
particular, in Euclidean space for (perturbed) gridded data, certain improved error 
estimates are already known to hold for functions within the native space itself (see, 

e.g., m). 

The goal in this paper is quite different from the desire to establish error estimates 
for functions possessing insufficient smoothness for admittance in the native space. 
In recent years, contributions in that direction has been provided by several authors, 
e.g., [T6l[T7, 12 for the sphere and 23, 3, 18J for the Euclidean case. 

2. The sphere. Let S'' = {a; G R'^'^'^ : \x\ — 1}. Then the geodesic distance 
between points x,y G S'' is d{x,y) = cos^^xy, where xy denotes the usual inner 
product of vectors in We let denote the normalised rotationally invariant 

measure on the sphere and define the inner product 

(/,5)l2(S<*) — / f{x)g{x)diy{x). 

Let II • ||L2(s<i) •= ( • J ■ )i,2^(s<i) ^^'^ -^2(8'^) denote the set of functions for which 
II " IIl2(S'') < ^'^^ be the polynomials of degree n in W^+^ restricted to the 
sphere, and let _ff„ = P„ n be the space of degree n spherical harmonics. Then, 
i2(§'*) has the decomposition 

n>0 

Let Y", . . . , Y2 be an orthonormal basis for iJ„. 

Related to S'' (we will see why shortly), we have the Gegenbauer polynomials 
C^^^(t) which are orthogonal on [—1, 1] with respect to the weight (1 — t^Y' ^1"^ . It is 
well known (Miiller for instance) that the following addition formula holds: 

C';^\xy)^Y.^^^x)Y^^y), 

with A = (i/2 — 1. The normalisation of the Gegenbauer polynomials is chosen so that 
there is no constant in the addition theorem. It is straightforward to see that Ci^' IS 
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the kernel of r„, the orthogonal projector from L2{§'') onto Thus, 



(T„/)(a;)= / /(y)C(^)(xy)dKy), 



for aU / e L2 (§''). 



The following lemma is a specialisation of a result in [TT] to the sphere. 
Lemma 2.1. For n> 0, 



for allfe L2{B'^). 



We will be considering interpolation using kernels of the form (j)(d{x,y)) where 
[0, tt] R. We will assume that the function (/) has an expansion 



(/.(d(x,y)) = ^a„C(^)(xy), 



n>0 



where a„ > 0, for n = 0, 1, 



and 



< 00. 



The first condition ensures that (j) is positive definite, and the second that it is con- 
tinuous. Our analysis will take place in the native space for </>, JC^, defined by 



:= \ f e US') : \\fU:= (j2a-'\\T„f\\ 



n>0 



< 00 



A pseudodifferential operator A on is an operator which acts via multiplication 
by a constant on each eigenspace iJ„: 



Ap„ = XnPn 



Pn e Hn, n^O,l,.. 



For more information on pseudodifferential operators on spheres see, e.g., [HI HO]- We 
call the sequence of numbers {A„},i>o the symbol of A. Let Sx denote the point 
evaluation functional at x, and, when it makes sense for the functional fj,, let Afi(f) = 
/z(A/). Let us denote by Y* the span of the point evaluation functional supported 
on Y. In Morton and Neamtu T¥] the authors give error estimates for the collocation 
solution of pseudodifferential equations on spheres. Here we attempt, initially, to find 
errors in pseudoderivatives of solutions to the interpolation problem. 

Proposition 2.2. Let he the (p-spline interpolant to f E 'K^ on the point set 
Y CZ . Let A be a pseudodifferential operator. Then, for each x E S'^, 

\Hf^Sl)ix)\< inf sup \Av{x)~f,{v)\\\f~SjU. 

p.eY' ^^^^ 

Proof Since /(y) — ^^(y) — 0, y eY, we have, for any coefficients Cy, y &Y, 



\A{f-Sl){x)\ = 



Aif-si)ix)^Y.^y(f(y)-sJ(y)) 



yeY 



A{f-Sj){x)-J2cyifiy)-Sj{y)) 



Wf-s] 



< sup 

lklL=i 



Av{x) - Cyv{y) 

V^Y 
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We now take the infimum over all functionals in Y* to obtain the result. □ 

In what follows we will need the pseudodifferential operator A to satisfy the 

following assumption: 

Assumption 2.3. For all n > 0, Xn — {n{d + n - 2))", for some s > 0. From 

Ditzian 6 , if A satisfies Assumption 12. 3i then for p G P„, 

for some E independent of n. 

From [ini Lemma 7] we have the following result. 

Lemma 2.4. Let Y be a finite set of points with fill-distance h{Y,§'^) < 1/{2N), 
for some fixed N G Z+ . Then, for any linear functional 7 on Pjv with 

sup \jp\ < 1, 

IIPlll-oo(S<i)=l 

there is a set of real numbers {by}yeY7 with "^y^yl^yl — ^' ■^^c/i that 
7P X! /o'' P& Pn- 



Now, for a fixed a; G S'^, let 



Then, 



7P=— — , for all pG Pat. 

tjAN 



sup II 11'^^' <1, 



so that, by the previous lemma, there is a set of coefficients {hy}y<=.Y, such that 



with X])^eyl^j/I — 2- Thus, with Cy = EX^by, for y £Y, we have 



where, 



Kp{x) ^^Cyp{y), forallpGPjv, (2.1) 

y(^Y 



^\cy\<2EXN. (2.2) 



We now arrive at the first main result of this section. 

Theorem 2.5. Let be the (j)-spline interpolant to f E CK^f,, on the point set Y C 
where h{Y,§'^) < 1/{2N), for some fixed N G Z+. Let A be a pseudodifferential 
operator with symbol {A„}„>o satisfying Assumvtion \2.3\ and 



dnXl^ttn < 00. 



?i>0 
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Then, for x G S'' 



|A(/ - Sj){x)\ <{1 + 2E)(Y, dnXlaA \\f - S, 



i>N 



'0 lU- 



Proof. Let is choose {cy}y^Y to be the coefficients described above. Let v G 
with \\v\\^ = 1. Then, 



inf |Aw(2:) — fJ,{v)\ < 



n>0 



f \nTnV{x) - ^ CyTnV{y) 
n>N ^ yeF 



by (123]). Thus, 



inf |Av(a;) < 



A„r„i;(a::) + ^ ^ CyTnv{y) 

n>N n>NyeY 

J2 + Z!l^aljll^"«llLoo(S'') 



n>N ^ j/ey ' 

using Lemma |2. II Hence, using (|2.2[) and the Cauchy-Schwarz inequahty. 



inf \Av{x) ~ fi{v)\ < 

tier' 



n>N 



n>N 



and the result foUows from Proposition 12.21 since ||w||0 = 1, and because {A„}„>o is 
an increasing sequence. □ 

Integrating the conclusion of the previous theorem over the sphere we easily obtain 

Corollary 2.6. Under the hypotheses of Theorem [KR 



l|A(/-5j)L,(S^) < {l + 2E)(Y,dnXlaA \\f - S, 



n>N 



by 



Before we give our improved error estimate we need to define a new space ^A<j> 



^A^:^ If e^K^: ||/||a0:= (^(A„a„)-2||T„/||2^(s,)) <^ 



Theorem 2.7. Let Sj be the (p-spline interpolant to f ^ "K^ on the point set Y C 
2)'^, where h{Y,2)'^) < 1/{2N), for some fixed N € 2+. Let A be a pseudodifferential 
operator with symbol {A„}„>o satisfying Assumvtion \2.3\ and 

n>0 
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Then, for f G !Ka0 and for all x G S"^, 

\f(x) - Sj{x)\ <{1 + 2E)(Y1 dnXW) ' ( dnan 

--^ AT / \ AT 



|A0. 



n>N ' ^n>N 

Proof. Firstly, using Proposition [TTl] and the Cauchy-Schwarz inequality, we have 

11/^ I I = (/ ^ T f)4> 

n>0 

< (E ^nWTnif - SJ)\\l^,,^ ' (E(^"«")-'rn/llL(S^) 

= ||A(/-5^)|U,(g.)||/||A0 
<{1 + 2E)(Y1 dnXW) ' 11/ - SjUWfU^, 
using Corollarv l2.6l Cancelling a factor of ||/ — 'SjlU from both sides yields 



n>N 



|A0. 



We can now employ the standard error estimate taken from Jetter, Stockier and 
Ward TU] (our Theorem 12.51 with A„ = 1 for all n) to give the required result. □ 

3. The Euclidean case. Our attention now turns to (/)-spline interpolants of 
the form 

where (f> : E_|_ — > R. We will conduct our analysis for positive definite basis functions 
4> € Li(R'^) whose Fourier transform satisfy, for some s > 0, 

+ i^r^^- < ^(x) < C72(i + (3.1) 

for some positive constants Ci and C2, for example, the Sobolev splines [7] or piecewise 
polynomial compactly supported radial functions of minimal degree |21j . The expo- 
sition contained in this section can be readily adapted to include the polyharmonic 
splines [5: as well. In that case, the (/)-spline interpolant must be augmented by a 
polynomial p with the extra degrees of freedom taken up by the side conditions 

where q is polynomial of the same degree (or less) as p. 

For a domain C M'' let L2(ri) denote the usual space of square- integrable 
functions on fl with inner product ( ■ , ■ )L2(n) ^-^d norm \\ ■ \\L2(n)- For k € Z+, the 
integer-order Sobolev space is defined by 

:= |/ e L2{R'^) : D"./ G i2(M'^) for all |a| < /c|, 
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with £)" understood in the distributional sense, which carries the inner product 

{L9)k ■■= {f,g)L2(R'^) + ifi9)k, 

where {f,g)k denotes the Sobolev semi-inner product 



if, a) 



k 



with associated semi-norm | ■ |fc :=(•,• )fc ■ "^^^ coefficients Ca have been chosen so 
that 



\a\=k 



- \x?K 



We can write the semi-norm, using the Fourier transform, in the alternative form 

\f\l = I |/(a:)P|:.p'=dx, 



which facilitates the definition of fractional-order Sobolev space, JCs, for s > 0, which 
has the semi-norm 



/ Ifixm^'dx. (3.2) 
The space !Ks is complete with respect to 

1 

niL(K-) + l/l')' ifsez+, 

1 

11/11 Lsj +1/1')' otherwise, 

and, whenever we have s > d/2, !Ks is continuously embedded in the continuous 
functions. The native space for </> satisfying is equivalent to 'Ks- 

We now wish to make local definitions of our function spaces, which we shall 
denote by 'Ksiyt). For s G Z+ the definition should be is clear. In what follows we 
also need the local fractional-order Sobolev spaces: 

:={/ e ^ys\m : \\f\\s.n := (ll/ll?.j ,o + l/ll^.) ' < ^ 

where |/|s,o is the local fractional-order Sobolev semi- norm obtained by rewriting (|3.2[) 
in an equivalent direct form, i.e., not defined through the Fourier transform of / (see, 
e.g., Adams [H p. 214]). For our analysis we find it more useful to exploit an equiva- 
lent wavelet representation for the local Sobolev norm [5j . 

To introduce this equivalent norm we stipulate that the bounded domain, 17, 
admits a local multiresolution of closed subspaces {Ki(f^)}n>o of L2(ri): 



n>0 

Cohen et al. fS] give sufficient conditions on fl to admit such a local multiresolution. In 
particular, for d — 2, those domains whose boundaries have certain piecewise Lipschitz 
smoothness are admissible. The following is an incidence of [5} Theorem 4.2]. 
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Theorem 3.1. Suppose fl C M is a bounded domain that admits a local multi- 
resolution {Vn{fl)}n>o for L2{^). For n > 0, let denote the orthogonal projection 
from L2{^) onto Wn{^) = Vn{^) V„_i(r2) with the convention that V-i{n) = {0}. 
For each s >0, let Ag be the pseudodifferential operator on Q. defined via 

n>0 

Then, there exists positive constants Ci and C2 such that, for all f G '}is{^), 
Ci\\Asf\\L,[n) < < C2||A,/|U,(o). 



Now, let us return to the task at hand. For (j) satisfying (|3.ip . we will denote the 
0-spUne interpolant on the point set Y by . The standard error estimate in this 
context is 

\f{x)-S^{x)\<Ch^-'/^f-Sj\U, (3.3) 

see [22] . If / is smoother (and satisfies some boundary conditions) we can get a better 
rate of convergence. 

We will exploit the fact that J{^(0), for < ^ < s, is an interpolation space 
lying between ^2(57) and J{s(ri) (see Bergh and Lofstrom [2l p. 131]). We can then 
use the standard interpolation theorem concerning the norms of operators bounded 
on the extreme spaces to infer a bound on the norm for the interpolation space. For 
further information on interpolation spaces the reader can consult, e.g., Bergh and 
Lofstrom We use the following interpolation theorem. 

Proposition 3.2. Let < fi < s. Further suppose that T : 'Hs{n) L2{fl), and 
T : !Ks(ri) — » (17) is a bounded operator. Then, 

l|r|L'H:,(0)^Mf,(0) < ll^llM^^f2)^L2(0)ll^llMj*(0)^M,(0)- 



Since we can write the JCs-norm in entirely direct form, we are at liberty to 
utilise Duchon's localisation technique [8] to enhance the standard error estimate (|3.3p . 
Therefore, if is bounded and satisfies an interior cone condition then, for / G !Ks(fi), 
s > d/2, and sufficiently small h = h{Y, ft), 

11/ - S^\\L2{n) < Ch^Wf - SjWs^n. 

Writing Tf — f — , we see, using the last proposition, that, for < < s, 

11/ - 5jiu,o < {ch^wf - sjUny-^'/'Wf - s^rj^ 

^Ch^-^^\\f-Sl\U,n. (3.4) 

We can now prove our main result of this section, which is a generalisation of 
that of Schaback [T5] . 

Theorem 3.3. Suppose C M"* is bounded, satisfies an interior cone condition 
and admits a local multiresolution. Let s > d/2 and let S^ be the (f)-spline interpolant 
to f £ J{s on the point set Y C fl. Suppose further that f G !K^, for s < v < 2s, and 
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that f is compactly supported in Q. Then there exists C > 0, independent of f and 
h — h{Y, fl), such that for all x ^ and sufficiently small h, 

\f{x)-Sl{x)\<Ch'^-'''^\\fU,. 



Proof. From Proposition 11.11 wc know that 

(/ ~ ^ Sj)s = 0, 

so that 

11/ - 5^11' = (/ - Sj,f)s < C{f - /),,n, 

where we have used the compact support of / in f2. Now, the equivalent norm from 
Theorem 13.11 gives us 

||/-5Xll?<C(A,(/-5X),A,/)i,(o) 

n>0 

and successive apphcations of the continuous and discrete Cauchy-Schwarz inequahty 
yields 



\\f-Sj\\l<Cj2\\2-^'^-'^^Q2if~Sl)h,^n)\\2-^Q2fh 



2(0) 



Tl>0 



= C\\A,,^M~Sj)\\L,^n)\\Kf\\L.in)- 

Thus, using the norm equivalence from Theorem 13.11 again together with 
have 

\\f-S^\\l<C\\f-S^hs-.,n\\f\l,n 
<Ch-'-^\\f~Sl\\ssAf\kn 
<Ch'^-'\\f-Sl\\s\\f\Un, 

and cancelling one power of ||/ — SY\\s gives 



\\f-Sj\\s<Ch-^-'\\f\\,,n. 

The result follows by substitution into the standard error estimate p.3p . □ 
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input into questions of pseudodifferential operators and alternative Sobolev space 
semi-norms. 

REFERENCES 

[1] R. A. Adams, Sobolev spaces, Academic Press, New York, 1975. Pure and Applied Mathematics, 
Vol. 65. 



10 R. A. BROWNLEE, E. H. GEORGOULIS, AND J. LEVESLEY 

[2] J. Bergh and J. LOFSTROM, Interpolation spaces. An introduction, Springer- Verlag, Berlin, 
1976. Grundlehren der Mathematischen Wissenschaften, No. 223. 

[3] R. A. Brownlee and W. Light, Approximation orders for interpolation by surface splines to 

rough functions, IMA J. Numcr. Anal., 24 (2004), pp. 179-192. 
[4] M. D. BuHMANN, Radial basis functions: theory and implementations, vol. 12 of Cambridge 

Monographs on Applied and Computational Mathematics, Cambridge University Press, 

Cambridge, 2003. 

[5] A. Cohen, W. Dahmen, and R. DeVore, Multiscale decompositions on hounded domains. 

Trans. Amcr. Math. Soc, 352 (2000), pp. 3651-3685. 
[6] Z. DiTZL'VN, Fractional derivatives and best approximation. Acta Math. Hungar., 81 (1998), 

pp. 323-348. 

[7] J. G. Dix AND R. D. Ogden, An interpolation scheme with radial basis in Sobolev spaces 
//=(R"), Rocky Mountain ,1. Math., 24 (1994), pp. 1319-1337. 

[8] J. DuCHON, Sur I'erreur d'interpolation des fonctions de plusieurs variables par les D"^-splines, 
RAIRO Anal. Numer., 12 (1978), pp. 325-334, vi. 

[9] W. FreeDen, T. Gervens, and M. Schreiner, Constructive approximation on the sphere. 
Numerical Mathematics and Scientific Computation with Applications to Geomathematics, 
Oxford Univ. Press, New York, NY, 1998. 
[10] K. Jetter, J. Stockler, AND J. D. Ward, Error estimates for scattered data interpolation 

on spheres. Math. Comp., 68 (1999), pp. 733-747. 
[11] J. Levesley, C. Odell, and D. L. Ragozin, Scattered data on homogeneous manifolds: the 
norming set approach, in Curve and surface fitting (Saint-Malo, 2002), J.-L. Merrien A. Co- 
hen and L. L. Schumaker, eds.. Mod. Methods Math., Nashboro Press, Brentwood, TN, 
2003, pp. 269-278. 

[12] J. Levesley and X. Sun, Approximation in rough native spaces by shifts of smooth kernels 

on spheres, J. Approx. Theory 133 (2005), pp. 269-283. 
[13] W. A. Light .vnd H. Wayne, On power functions and error estimates for radial basis function 

interpolation, J. Approx. Theory, 92 (1998), pp. 245-266. 
[14] T. M. Morton and M. Neamtu, Error bounds for solving pseudodifferential equations on 

spheres by collocation with zonal kernels, J. Approx. Theory, 114 (2002), pp. 242-268. 
[15] C. Mi'iLLER, Spherical harmonics, vol. 17 of Lecture Notes in Mathematics, Springer- Verlag, 

Berlin, 1966. 

[16] F. J. Narcowich, R. Schaback, and J. D. Ward, Approximations in Sobolev spaces by kernel 

expansions, J. Approx. Theory, 114 (2002), pp. 70-83. 
[17] F. J. Narcowich and J. D. Ward, Scattered data interpolation on spheres: error estimates 

and locally supported basis functions, SIAM J. Math. Anal., 33 (2002), pp. 1393-1410 

(electronic). 

[18] , Scattered-data interpolation onR"; error estimates for radial basis and band-limited 

functions, SIAM J. Math. Anal., 36 (2004), pp. 284-300 (electronic). 
[19] R. Schaback, Improved error bounds for scattered data interpolation by radial basis functions. 

Math. Comp., 68 (1999), pp. 201-216. 
[20] S. L. Svensson, Pseudodifferential operators - a new approach to the boundary problems of 

physical geodesy, Manuscripta Geod., 8 (1983), pp. 1—40. 
[21] H. Wendland, Piecewise polynomial, positive definite and compactly supported radial functions 

of rninirnal degree. Adv. Comput. Math., 4 (1995), pp. 389-396. 
[22] Z. M. Wu ,'\nd R. Schaback, Local error estimates for radial basis function interpolation of 

scattered data, IMA J. Numer. Anal., 13 (1993), pp. 13-27. 
[23] J. YoON, Interpolation by radial basis functions on Sobolev space, J. Approx. Theory, 112 

(2001), pp. 1-15. 



